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Abstract 

In this paper the modification of the method conventionally used for 
the modeling of the massive fermions production and decays is pro- 
posed. The step by step algorithm is presented. Under the strict 
conditions the proposed method of modeling allow distinctly raise the 
efficiency of the computations. 
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Simulation of massive fermion production (let us denote it F) with sub- 
sequent decay frequently occurs in High Energy Physics. Indeed, the angular 
distributions of the fermion decay products are very sensitive to its polariza- 
tion. 

The most obvious way is the calculation of the squared amplitude of the 
whole process including the fermion decays. However such way may be ac- 
companied by some difficulties. First, the calculation of the amplitude of the 
whole process (including the decay chain) could be rather cumbersome. More 
over the corresponding final state can appear without the massive fermion 
contribution. In other words the final state of the process could be the result 
of the additional diagrams not containing F at all (see Fig. [1]). 
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Figure 1: The example of the diagrams describing the process with the same final 
state. Diagram a) corresponds to the massive fermion F resonance contribution. 
Here F decays into bW . Diagram b) has the same final state, but not contains 
fermion F. 

However, one can assume that the invariant mass of the F decay products 
is close to fermion mass mp- In this case the contribution of non-resonance 
diagrams (see FigHb) will be notably suppressed. 

Fairly often the problem could be simplified with the narrow resonance 
approximation, i.e. the intermediate fermion is considered to be «on the 
shell» {pj^ = rri^p^ pp is 4-momentum of the fermion F). Such approach 
allows to simplify the calculations and split the simulation into two steps: 

1. the production of the massive polarized fermion; 

2. the subsequent decay of the polarized fermion 

Thereby the squared amplitude of the process can be written in the «factor- 
ized» form: 

|^|2 ^ |^Prod^^^)|2 . lA'^'^s)]^ (g) $(S). (1) 
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Here A^^°'^{s) and A^'^'^{s) are the amplitudes describing the production and 
decay of the polarized fermion, s denotes the spin of the fermion {{spp) = 0) 
and is the factor of «spin transfer* from production to decay. 

The problem posed in such way used to be solved within the method 
of spiral amplitudes In this method the amplitudes of the fermion F 
production and decay calculated depending on the different helicities (Aj). 
Than the total amplitude of the process can be performed as the following 
sum: 

\A\^ = J2 A^™''(Ai)A^'='=(Ai) ■ (A^™''(A,))+(A^'='=(A,))+, (2) 

The equivalent approach was proposed by Jadach and Was [2]. They 
showed that expression ([2]) can be rewritten as follows: 

\A\^ = |yl^™'^p ■ \Aj^'fil + HV). (3) 

Here A^'^"'^ and Aq^'^ are the amplitudes of unpolarized fermion production 
and decay. H and V are so-called polarization vectors. They are determined 
in the fermion rest frame and contain information about fermion spin. 

|^Prod^^^|2 ^ (1 + (i/s)) (4) 

= K^f + (5) 

One example is given below. Accordingly to ((51) the matrix element 
squared for the t quark decay to hW^ with the following decay of boson 
to can be written as follows: 



I M |2 = iPbPu){PtPi) 

' ' iPl-M'^y + Tl.M'^ 



(PtPe) 



(6) 

'ly L yPtPi) J 

where pt,Pb,Pi,PnU are the momenta of the t and b quarks and final leptons, 
Tw and are the total decay width and the mass of the VF-boson. Hence 
one can derive the polarization vector for this decay 

As an example for expression ((31) the process ud tb, t —>■ biu is 

considered (see Figl2l). 

Using (JG]) the explicit expression for the matrix element squared of this 
process can be written in the form (compare to ([3])) as follows: 

4 {PuPb){PtPd) 



\M\' oc [g 
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Figure 2: Feynman diagram for the process ud bivb 

where and are directions of and (i-quark momenta in t-quark rest 
frame. 

Fairly often the expression ([3]) is preferable in the numerical calculations 
since it reduces the quantity of logical operation. Such conventional method 
of simulation uses discarding technique (the reject-and accept method) [3] 
and supposes realization of the following algorithm (see details in [2]): 

1) to calculate the filal momenta of the particles from the fermion F pro- 
duction process; 

2) according to a given kinematics to evaluate the polarization vector H 
in F rest frame; 

3) to perform simulation of the F-fermion decay and to fix the decay 
products kinematics; 

4) to evaluate the decay polarization vector V (in F rest frame); 

5) to calculate the additional weight W = (l + HV) and using a discarding 
technique to reject or accept an event; 

6) if the event is discarded than return to step 3). 

However using this algorithm one can expect the increase of computation time 
in case of complicated expressions for the massive fermion matrix elements. 

In this Note the modification of this algorithm is proposed. Very often 
such modification allows to raise an efficiency of the numerical calculations. 
The basic condition of the modified algorithm application is \V\ = const, 
i.e. the vector V absolute value must be independent of the fermion decay 
simulation results. 

One should notice that the value of the additional weight W = {1 + HV) 
(see step 5 of the algorithm) depends on the direction of V only. Besides the 
different kinematics of the decay can lead to the same polarization vector V . 
This ideas are the basic for our modified algorithm. 

Thus we offer to chose correctly the polarization vector V (actually the 
direction of V) before the fermion decay simulation. Correctly means in 
order to accept an event. Hence the algorithms steps starting from step 3) 
modifies: 



3 



3as) chose the V direction using an discarding technique with the weight 

W = l + HV; 
4as) to simulate the massive fermion decay ; 

5as) to evaluate the polarization vector v corresponding to the F decay 
kinematics (mote, that in general v ^ V since the vectors may have 
different directions); 

Gas) to rotate the reference system so that vector v coincides with V in F 
rest frame . 

Let us point out that for each accepted event the simulation of fermion 
decay made just once. This fact is for sure an advantage of the proposed 
algorithm. Note, that the method described above can be laso used in case of 
two fermions production with a subsequent decays. In this case the additional 
weight takes the form: 

W2 = l + hivi + h2V2 + hikv\v^. (9) 

where (like in ([3])) two vectors hi^2 and tensor hik are determined by produc- 
tion kinematics. The simulation of such two fermions decays and the turning 
of the reference systems should be made independently. 
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